Evolutional Entanglement in Nonequilibrium 

Processes 

V.I. Yukalov 

Bogoluhov Laboratory of Theoretical Physics 
Joint Institute for Nuclear Research, Duhna 141980, Russia 



Abstract 

Entanglement in nonequilibrium systems is considered. A general defini- 
tion for entanglement measure is introduced, which can be applied for charac- 
terizing the level of entanglement produced by arbitrary operators. Applying 
this definition to reduced density matrices makes it possible to measure the 
entanglement in nonequilibrium as well as in equilibrium statistical systems. 
An example of a multimode Bose-Einstein condensate is discussed. 

1 Introduction 

The concept of entanglement [1] is believed to play an important role in quantum 
information processing and quantum computing, because of which much efforts has 
been devoted to quantifying entanglement [2-8]. The latter is usually measured by 
a kind of reduced or relative entropy and is considered for bipartite systems [9,10]. 
The choice of entropy as a measure of entanglement looks reasonable since entropy 
characterizes the complexity of both equilibrium and nonequilibrium systems [11]. 
However, measuring entanglement by means of a reduced or relative entropy is 
relatively straightforward only for rather simple bipartite systems. Moreover, there 
does not exist a measure of many-body entanglement for multipartite systems. It is 
even less clear how to quantify entanglement in nonequilibrium statistical systems. 

The aim of the present communication is to introduce a general entanglement 
measure that would be valid in arbitrary cases and to illustrate its application to 
a nonequilibrium system. As an example of the latter, a multimode Bose-Einstein 
condensate of trapped atoms is considered. 



2 Entanglement Measure 

Before going to physical applications, it is useful to define what actually entangle- 
ment means from the general mathematical point of view. As far as one usually 
speaks about the entanglement of some states, it is necessary, first, to specify what 
states are to be entangled. 

Consider a set of objects, called parts, which are enumerated by an index i = 
1,2, ... ,p. Each part is characterized by a Hilbert space of single-partite quantum 
states 

n, = C{\n,>}, (1) 
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being a closed linear envelope of a single-partite basis { \ni >} of quantum states 
\ni>. The space of composite-system quantum states is a subspace 



(2) 



of the p-fold tensor product Ti^ of spaces (1). In particular, Ti may coincide with 
7i^. But 7i becomes a subspace of if the tensor-product space is complimented 
with some selection rules restricting admissible states. For example, such selection 
rules may consist of the requirement of specific symmetry properties. A subspace of 
a tensor- product space can be called incomplete tensor-product space [12]. 

From the composite-system space 7i, we separate out the set T) G Ti. oi disen- 
tangled states, thus, obtaining the disentangled set 



The states of V have the structure of tensor products of (fi G Hi, but V does 
not include linear combinations of such products. This principally distinguishes 
the disentangled set T> from the total space H. The latter, in addition to the 
tensor products of quantum states, contains as well their linear combinations. The 
compliment 7i \ P composes the set of entangled states. 

The meaning of the word entanglement is that the disentangled states from 
V are, by means of some transformation, converted into the entangled states of 
'H\V. Transformations are accomplished by operators. Hence, mathematically, 
entanglement implies that an operator A acting on T) transforms it into 'H\T>. 
Thus, the concept of entanglement, to be mathematically correct, must include 
the definition of disentangled states, to be entangled, and the specification of an 
operator, producing this entanglement. 

Let a bounded operator A be defined on Ti,. Because it is bounded, it possesses 
a finite norm 



where | l^^l |>t is a vector norm of </? on H.. In addition to the norm A on H., we may 
define the norm of A on D, that is. 



It looks that to quantity the entanglement caused by an operator, we shall need 
to invoke the operator norms. This understanding comes from the fact that the 
operator norms are involved in measuring the order associated with operators, which 
is characterized by the operator order indices [13]. These indices can be introduced 
for arbitrary operators. In particular, they can be defined for density matrices, thus, 
describing ordering in physical systems [14]. 

The entanglement, produced by an operator, should be described by comparing 
the performance of this operator on the disentangled set with the action of a related 
operator that does not entangle the states of V. A nonentangling operator should 




(3) 
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sup 

ll¥'llw=l 



A\\v = sup 

Il/lb=l 
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have the structure of the tensor product of operators acting on Hi- To this end, we 
construct the product operator 

in which 

A\ = const Tr{^^._^.} A 
and the form (4) is chosen so that to satisfy the normahzation condition 

Tt-hA^ Ttv A^ . 

Comparing the action on V of an operator A with that of its nonentanghng coun- 
terpart (4), we define the entanglement measure 



\A\ 



<A) ^ log iri^ > (5) 
Id 



where the logarithm can be taken with respect to any base, say, to the base 2. 

It is straightforward to show that the entanglement measure (5) possesses the 
properties that are natural for being such a measure. First of all, this measure is 
semipositive, 

e{A) > . (6) 

Second, the measure is continuous in the sense that if for any operator AonH there 
exists a family {^(i)} of operators A{t), parametrized with i e R, so that 

\\A{t)\\v ^ \\A\\^ {t^O), 

then 

s{A{t))^eiA) (t^O). (7) 

Third, a nonentangling operator, having the structure of a tensor product A®, does 
not produce entanglement, 

e (a®) = . (8) 



The latter property may ne generalized to the case when A = (BuPi^Af is a linear 
combination of the operators Af such that 



V — ^ V 



In that case, one has 

e [®u PuAf) = . 

Fourth, the entanglement is additive, which means the following. Let A = ^,^A^, 
then 

e = J2 ^(^v) . (9) 
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Finally, the measure is invariant under local unitary operations Ui, such that U^Ui — 
1, for which we have 

s{®UUU®UU^)=eiA). (10) 

In this way, the entanglement measure (5) can be defined for an arbitrary op- 
erator. In physical applications, one may consider the entanglement produced by 
any operator from the algebra of observables, for instance by a Hamiltonian, by a 
spin operator, and so on. One may also investigate the entanglement caused by 
statistical operators and by density matrices. 

As an illustration, wc may consider a simple case of a bipartite system, when 
Ti. = Til ® 7^2- Let each Tii be a separable Hilbcrt space of dimension di, which is 
a span of an orthonormal basis { \ip'^^ >}. Consider entanglement realized by a von 
Neumann statistical operator 

Pbp=\BP><BP\, \BP>en, 

describing a pure statistical state of the bipartite system. The wave function of a bi- 
partite system can be presented (see e.g. [14]), involving the Schmidt decomposition, 
as the biorthogonal sum 

d 

\BP> = Y1 C„ \iPn > ® > , 

n=l 

in which d = min^ di and Y.t=i — 1- Following the procedure, described above, 
we have 

d 

P\ = ^H,^i PBP = II |Cn 
n=l 

The corresponding norms are 

WpbpWv ^WPiWui =sup|c„|^. 

n 

Then for the entanglement measure (5), we get 

£{pBp) = -logsup|c„|^ . 

n 

This varies in the interval 

< e {pBp) <^ogd . 

Maximal entanglement occurs when |c„p = \/d. One usually quantifies entangle- 
ment in a pure bipartite system by the reduced von Neumann entropy 

d 

Si^ = -Trw, p\ logpi = -Yl \Cn\^ log \Cn\^ . 

n=l 

The latter coincides with measure (5) for the maximally entangled state, when 
S}^ — s{pbp) = logo?. Note that for a mixed bipartite state the reduced entropy Slf 
is, generally, different for i = 1 and i = 2, thus, becoming not well defined as an 
entanglement measure, while e{p) is always well defined. 
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3 Evolutional Entanglement 



In general, if the considered operator A = A{t) depends on time t, the entanglement 
measure (5) will be a function of time e{A{t)), displaying the temporal evolution of 
entanglement. To study the evolutional entanglement for physical systems, we may 
employ the reduced density matrices. 

Consider a set x'^ = {xi, X2, ■ ■ ■ , Xp} of variables characterizing a physical system. 
Each Xi G Xi pertains to a characteristic space Xi. A p-order reduced density matrix 

= (11) 

is a matrix with respect to the variables and x^, The matrix elements being 

Pp {xP,xP,t) = Tr^ iPixi) . ..i/;{xp)p{t)i^\xP) . . . V'^(xi) , (12) 

where the trace is over the Fock space, i/j{x) is a field operator, and p{t) is a statistical 
operator. The first-order density matrix is 

p\{t) = [p,{x^,x^,t)] , (13) 

with the elements 

pi {x,x,t) = Tijr ■ip{x)p{t)-ip''{x) = < ^|J\x) ^jJ{x) > . (14) 

One often calls the matrices (11) and (13) the p-particle and single-particle density 
matrices. This in no way means that some concrete particles are compulsory sep- 
arated out of the system, but just designates the order of the matrices. We may 
keep in mind a system of N indistinguishable particles, with the same character- 
istic space X. Then a p-order density matrix describes correlations between any p 
particles from the ensemble of N identical particles. 
The trace of p\{t) over the single-partite space (1) is 

Tr^ii P\it) = J2 IPiWI n'i> = I pi{xi,Xi,t) dxi^N . (15) 

Similarly, the trace of Pp{t) over Ti is 

Trn Ppit) =JPp (x^, x^, t) dx^ = . (16) 

The relation between the single-partite and p-partite density matrices reads 

The product operator (4) becomes 

= (TV - rt! TV. «<- pW) ■ (18) 
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Calculating the norms of density matrices, over the disentangled set (3), we may 
define the single-partite basis >} as formed by the eigenvectors of p\{t). Then 
the latter can be presented as the diagonal expansion 

p\{t) = Y.DnAt)\ni><ni\. (19) 

rii 

For the entanglement measure (5), we have 

eip it)) = log (N-Py-N''\\PMh .20) 

This describes the level of entanglement between any p parts from the system of N 
parts. 



4 Multimode States 

To illustrate the calculation of the entanglement measure (20), let us consider the 
case of a multimode coherent system whose p-partite density matrix can be reduced 
to the form 

Pp''^'> = (N-p)\ ^ \n...n><n...n\ , (21) 
in which the fractional mode populations Wn{t) satisfy the properties 

0<Wn(t)<l, Y^Wn{t)^l. (22) 

n 

As a physical application, we may keep in mind a system of A^ trapped atoms in Bose- 
Einstein condensate (see reviews [15-17]). At low temperature, when practically all 
A" atoms are condensed, the system is in a coherent state [17,18]. Several modes 
in such a system can be created in different ways, e.g., by localizing atomic clouds 
in a multiwell potential formed by an optical lattice [19,20], by mixing falling and 
reflected atomic wave packets [21], and by other means. A controllable method 
of exciting topological coherent modes in a trapped Bose-Einstein condensate is 
by modulating the trapping potential with the help of resonant alternating fields 
[22-24]. 

For the density matrix (21), we have 

A^! 

"^^^^^"^ " iN-p)\ ^^P^"^^) ' 
and for the product matrix (18), we find 

A^' 

IW)II^=^]^3^ sup<(t). (24) 

Therefore, the entanglement measure (20) becomes 

^Mt)) = (1 -P) log sup Wn{t) . (25) 
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If the number of modes is m = J2n 1) then the maximal entanglement happens when 
all Wn{t) are equal with each other, that is, are equal to 1/m. The entanglement 
measure (25) can take the values in the interval 

< e{pp{t)) < (p - 1) log m . (26) 

The temporal behaviour of measure (25) is prescribed by the evolution of the frac- 
tional mode populations Wn{t). 

In order to be absolutely concrete, let us present an example explicitly demon- 
strating the evolution equations for Wn{t). For this purpose, let us consider the 
resonant generation of topological coherent modes in a trapped Bosc-Einstein con- 
densate [22-24]. An ensemble of coherent atoms is described by the Gross-Pitaevskii 
equation [16-18]. The stationary solutions to the latter define the topological coher- 
ent modes with a discrete spectrum of energies En- The condensate is subject to the 
action of external alternating fields with the frequencies tuned close to some of the 
transition frequencies Umn = {Em — En)/h. Suppose, for concrctcness, that there are 
two resonant fields whose frequencies are tuned to the transition frequencies UJ21 and 
U!s2, with the detunings A21 and A32, respectively. In the resonance approximation, 
the Gross-Pitaevskii equation can be reduced [22-24] to the system of equations for 
the mode amplitudes c„(i) which define the fractional mode populations 

Wn{t) = |c„|' . (27) 

In the considered case of two resonant fields, connecting three coherent modes, we 
have three mode populations, wi{t), W2{t), and W3{t). The evolution equations for 
the latter can be obtained from the equations for the mode amplitudes c„(i). The 
equations for the mode populations involve the ladder variables 

hi = 2cl 02 exp{i(A2it + 712)} , /12 = C3 exp{i(A32t + 723)} , 

hs = 2c* C3 exp{i(A32 + A2i)t + ^(712 + 723)} , (28) 

in which 7^„ are initial phases. For the mode populations (27), we derive the 
evolution equations 

^ = -bi2{h,-h) , 

^ = ^^23 - /12) - I bu {hi - hi) , 

^ = - ^ ^23 {h; - h2) , (29) 

where b„in are the transition amplitudes due to external alternating resonant fields. 
And the ladder variables (28) satisfy the equations 

. dh2_ ^ [ai2W2 - a2iwi + (ai3 - "23) + A21] - &12 {^2 -^"i) + 7; ^23 h , 
dt 2 

i = -/i2 [oL2zWz - a32W2 + {a2i - 0:31) Wi + A32] - 623 ('"^3 - '"^2) - 2 h , 
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i = -h^ [ai2,W2,-aziWi + {ai2-a2,2) u'2+A32+A2i] - - 612 ^2 + - ^23 hi , (30) 

where a^n are the transition amphtudes corresponding to the interatomic interac- 
tions. The system of equations (29) and (30) contains nine real-valued equations. 
However, by definitions (27) and (28), not all variables Wn{t) and hn{t) are indepen- 
dent. Thus, we have the following relations 

= 4-^1 ^2 , |/l2p = 4^2 , 1/131^ = 4^3^1, 

hi /l2 = 2^2 /l3 , Wi^-W2 + W^^l. 

As a result, Eqs. (29) and (30) can be reduced to an effective four-dimensional 
dynamical system. Solving these equations, we can find the temporal behaviour of 
the entanglement measure (25). 

The solution of the system of nonlinear equations (29) and (30) requires numeri- 
cal calculations. It is not trivial even for a simpler two- mode case [22-24], displaying 
such interesting effects as mode locking [22,24] and critical dynamics [23-25]. The 
situation is essentially more complicated for the three-mode case, considered here. 
In the frame of a brief communication it is impossible to discuss the details of the 
temporal properties of the mode populations (27), which will be done in a separate 
publication. We may only mention that the behaviour of Wnit) is relatively simple 
when \hmn\ <S 1, being a type of Rabi oscillations. Respectively, the entanglement 
measure (25) will oscillate with time. But with increasing |&rnn|, the temporal evo- 
lution of this measure becomes much more intricate. 



5 Conclusion 

A general definition for entanglement measure is advanced, which is valid for quan- 
tifying entanglement realized by an arbitrary operator. We concentrate here on 
discussing the case of nonequilibrium statistical systems, for which it is reasonable 
to study the entanglement caused by reduced density matrices. As an illustration, 
it is explicitly shown how the entanglement measure can be calculated for the case 
of a resonant Bose-Einstein condensate with several generated topological coherent 
modes. The resonant generation of these modes and their temporal behaviour can 
be governed by external modulating fields. Consequently, the evolutional features 
of entanglement can also be regulated. This suggests new possibilities for quantum 
information processing and quantum computing. 
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